The Arches Cluster - Evidence for a Truncated Mass Function? 
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ABSTRACT 

We have analyzed high-resolution, adaptive optics (AO) HK observations of 
the Arches cluster obtained with NAOS/CONICA. With a spatial resolution of 



84 mas, the cluster center is uniquely resolved. From these data the present-day 
mass function (MF) of Arches is derived down to ~ 4M . The integrated MF as 
well as the core and 2 nd annulus MFs are consistent with a turn-over at 6 — 7 M & . 
This turn-over indicates severe depletion of intermediate and low-mass stars in 
the Arches cluster, possibly caused by its evolution in the Galactic Center envi- 
^ \ ronment. The Arches MF represents the first resolved observation of a starburst 

in ' 

tions for stellar population synthesis modelling of extragalactic starbursts, the 
derivation of integrated properties such as the total mass of star clusters in dense 

Q-f chemical enrichment during starburst phases. 
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1. Introduction 



cluster exhibiting a low-mass truncated MF. This finding has severe implica- 



environments, the survival of low-mass remnants from starburst populations, and 
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The Arches Cluster is one of the rare compact starburst clusters forming in the Milky 
Way today. At a projected distance of only 25 pc from the Galactic Center (GC) it forms 
a unique template for star and cluster formation in the nuclei of external galaxies. Even 
at an age of only 2 Myr, its elongated shape evidences the strong tidal forces acting to 
disrupt the compact starburst cluster despite its large central density of 3 x 10 5 M pc~ 3 
(Figer et al. 1999). N-body simulations suggest that clusters in the GC environment will 
dissolve on timescales as short as 10 Myr (Kim et al. 2000; Portegies Zwart et al. 2002). 
This leaves us with very few clusters available to study nearby spatially resolved starburst 
clusters, especially at the Milky Way's present low cluster formation rate. 

The Arches cluster is known to be heavily mass segregated, with a flat present-day MF 
(r ~ where the Salpeter slope is V = —1.35) in the cluster core (r core = 0.2 pc at a distance 
of 8 kpc) (Stolte et al. 2002; Figer et al. 1999). The concentration of high-mass stars in the 
cluster core causes the integrated MF to be flattened with respect to a Salpeter MF with 
a slope of T ~ —0.9 to — 1.1 depending on the radial selection of the cluster area (Stolte 
2003). In this letter, we derive the present-day MF of the Arches cluster from AO data with 
a spatial resolution of 0. 084. This comprises the currently deepest photometry of the cluster 
core with the highest spatial resolution, which allows us to derive the present-day MF down 
to ~ 5 M Q in the cluster core (r < 0.2 pc) and 2 M at larger radii. 



2. Observations and data reduction 

NAOS-CONICA (NACO) is the ESO VLT facility AO system with its infrared camera. 
NAOS is a 189 element Shack-Hartmann AO system with two wavefront sensors operating 
in the visual and near-infrared (Rousset et al. 2003). CONICA is a 1 to 5 micron imaging 
camera and spectrograph designed for diffraction limited observations with the VLT 8m 
mirror (Lenzen et al. 2003). 

The Arches cluster {a = 17M5m50s, 5 = -28°49 28 (J2000)) was observed during 
the NACO commissioning phase in March 2002 to probe the AO crowded field performance. 
The medium resolution camera (O'.'027/pixel) optimised for if-band observations was used, 
covering a field of view of 27". A foreground star with V ~ 16 mag close to the cluster 
served as the guide source. Moderate Strehl ratios of 14 % in H-h&nd and 20 % in if-band 
were achieved, resulting in a resolution of 0'.'084 in both filters. Integration times were 14 
min in H and 7 min in K s . The detection limits on the combined images shown in Fig. 1 
were K s < 20.5 and H < 22 mag. 

Standard data reduction was performed under IRAF, images were combined using driz- 
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zle with 2x2 oversampling to enhance source detection, and PSF fitting photometry was 
obtained with DAOPHOT. The isoplanatic patch was large during the observations, such 
that no elongation of the PSF at larger distances from the guide star due to anisoplanatism 
is observed. HST/NICMOS observations of the Arches cluster (Figer et al. 1999) were used 
to probe the stability of the photometric correction and derive a uniform photometric cali- 
bration over the NACO field. The calibration procedure is described in detail with respect 
to the Gemini/Hokupa'a data set in Stolte et al. 2002. In the case of the NACO data, the 
zeropoint was constant and no color terms were found. In the following, all magnitudes are 
shown in the HST/NICMOS system, with K corresponding to the m205 NICMOS broadband 
filter, and H to ml60. All isochrones are transformed into the NICMOS system accordingly 
(Brandner et al. 2001). 



2.1. Artificial star simulations 



The completeness in the crowding limited cluster field was derived from the recovery 
fraction of artificial stars. In order to avoid changes in the local stellar density, only 200 stars 
were inserted per frame on a total of 100 frames, or 20000 stars. The same stars randomly 
inserted into the i^-band frame were transformed to instrumental iJ-band magnitudes using 
the characteristic Arches main sequence color of H — K = 1.63 mag. This is important 
because only stars detected in both H and K simultaneously enter the MF calculation. PSF 
photometry was derived in the same way as the original photometry. Recovered stars were 
selected by uncertainty and color in the same way as real stars entering the MF. The recovery 
fraction is determined for each annulus individually, such that increasing crowding effects at 
smaller radii are taken into account. In the cluster core, the stellar density is significantly 
higher than in the outer annuli. The simulations showed that inserted bright stars influence 
the recovery fraction of fainter stars severely. Therefore, the core field with the highest stellar 
density was excerpted from the NACO field. On this excerpt, only 10 stars within bins of 1 
mag were added on 100 artificial frames, leading to 1000 stars per magnitude bin. For the 
core number counts (r < 0.2 pc), only the recovery fraction from the core simulation was 
used. 



3. The Arches present-day MF 

Prior to the mass determination, stellar magnitudes and colors were corrected for the 
systematic variation in extinction observed over the Arches field (Stolte et al. 2002). Fore- 
ground and background sources were rejected by applying a color cut as shown in Fig. 2. 
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The increasing photometric uncertainty towards fainter stars as indicated by artificial star 
tests was taken into account by widening the color selection towards fainter magnitudes. 
Due to the location in the dense GC region, field contamination increases towards fainter 
magnitudes. This is apparent in the increase of blue foreground and red background sources 
in the CMD below K = 15.5 mag (14 M ), while contamination at higher masses is negli- 
gible. A dense concentration of bulge stars can be seen at the faint end, (K > 18 mag), 
which is beyond the 50% completeness limit in the MF and does not affect the MF fit. Field 
contamination becomes particularly important beyond a radius of 10". As no background 
field is available, the MF is limited to radii r < 10" or 0.4 pc, twice the core radius of 0.2 
pc (Stolte et al. 2002) to minimize field contamination. As some field star component is still 
present in the radial and color selection, this results in a conservative upper limit to the 
flattened MF. 

The present-day MF was derived from the corrected CMD (Fig. 2) converting if-band 
magnitudes into masses using a 2 Myr Geneva main sequence isochrone (Lejeune & Schaerer 2001) 
with solar metallicity 5 . A distance modulus of 14.52 mag and foreground extinction of 
Ay = 25.2 mag were applied to the isochrone, and reddened luminosities were transformed 
into the NICMOS system using color equations derived from comparison of five NICMOS 
standards and reddened stars in the Arches field with 2MASS. Present-day masses were used, 
as present-day properties are calibrated against observations while the re-calculation to ini- 
tial masses is heavily dependent on the mass-loss model. No correction was made for stellar 
evolution, which affects only stars above 50 M (log M/M Q > 1.7), to avoid uncertainties in 
stellar evolution modelling of very massive stars. 

Stars were counted in magnitude intervalls of A log M/M Q = 0.2. In order to minimize 
binning effects, ten MFs were calculated with bins shifted by 5\ogm = 0.02. The points 
from all ten calculations are included to yield the present-day MF presented in Fig. 3. 
The advantage of this method lies in the fact that features in the MF are not erased as 
a consequence of binning, and that the resultant MF fit is statistically more significant. A 
linear least squares fit was performed to yield the slopes. The true uncertainty in the derived 
MF slope is larger than the formal fit uncertainty, and is dominated by the choice of the 
isochrone, the photometric uncertainties of individual stars, and individual extinction. We 
estimate the true uncertainty to T ± 0.15 from experimenting with different fitting ranges 
and isochrone parameters (metallicity, age; see Stolte 2003 for details). 



5 Thc metallicity in the high-mass component of Arches was recently found to be solar(Najarro et al. 2004), 
while supersolar metallicity isochrones were used in earlier studies, see discussion in Figer et al. 1999, Stolte 
et al. 2002 and Stolte 2003. 
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The most striking feature in the present-day MF of Arches (Fig. 3) is a change in the 
sign of the MF slope at ~ 6M indicative of a turn-over (TO). The MF has therefore been 
fitted in the mass range 6 < M < 65 M (log M/M© > 0.8). The MF is with a slope of 
T = —0.86 moderately flattened with respect to a Salpeter value of T = —1.35. This slope is 
in very good agreement with previous values of T = —0.77 derived from HST/NICMOS data 
(Figer et al. 1999; Stolte et al. 2002), and T = —0.82 derived from Gemini/Hokupa'a data 
(Stolte 2003). Below M < 6 M (log M/M Q < 0.8) the MF deviates from a power-law, but 
appears to decrease towards lower masses between 0.6 < logM/M & < 0.8. Around 4M , 
the increase in the apparent MF is likely caused by the increasing contribution of field stars. 

As field contamination increases with decreasing mass especially for masses below ~ 
10 M Q (K = 16.3 mag), the observed TO will be enhanced when field stars are subtracted. 
As an upper limit to the field star contribution we consider stars close to the edge of the 
NACO field, at cluster radii r > 13" (0.5 pc). The corresponding apparent MF, scaled to 
the same area (Fig. 3, dash-dotted line), suggests that the population is dominated by field 
stars or segregated low-mass cluster members below 3M , while the high-mass end above 
15 Mq is entirely dominated by cluster stars. Between 3 and 15 M , the range where the 
TO occurs, field stars may contribute as much as 50%. The discrepancy between the high- 
and the intermediate mass regime in the Arches MF will therefore become more pronounced 
as soon as proper field star estimates are available. The TO in the MF is supported by the 
fact that no evidence for such a feature is observed in the outer field population, as would 
be expected if field star contamination introduces this feature into the MF. 

Two possibilities can be envisioned to explain the flattening in the Arches MF. The 
MF can either be truncated at the low-mass end, which implies that the dense GC envi- 
ronment has severe impacts on the local initial MF (IMF), or the tidal forces acting on 
the cluster caused rapid dynamical segregation on timescales as short as 2 Myr. Such dy- 
namical evolution would transport massive stars into the core and low-mass stars to larger 
radii. The dynamical timescales for an Arches-like cluster are suggested to be as short as a 
few Myr from N-body simulations (Kim et al. 2000; Portegies Zwart et al. 2002), such that 
dynamical disruption has likely contributed to the observed spatial distribution of cluster 
stars. 



4. The MF turn-over in two annuli 

In the case of dynamical segregation, the apparent high-mass TO would be a feature 
of the cluster core, where the massive stars accumulate, while stars with lower masses are 
ejected. The median mass is expected to decrease with distance from the cluster center. If 
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the median mass is evidenced in the TO, the TO mass should also shift towards lower masses 
at larger radii. The core (r < 0.2 pc) and 2 nd annulus (0.2 < r < 0.4 pc) MFs are compared 
in Fig. 4. While the core and 2 nd annulus MFs significantly vary in shape, the TO is present 
at logm ~ 0.84 or 7M Q in the core and logm ~ 0.8 or 6 M in the 2 nd annulus, in both cases 
above the 50% completeness threshold, i.e. only marginally higher in the core. Although the 
TO appears much more pronounced in the core, the TO mass is not shifted significantly 
to lower masses in the 2 nd annulus, as would be expected from dynamical evolution alone. 
The constant TO mass in the inner and outer cluster supports the suggestion of a heavily 
low-mass depleted IMF in the Arches cluster. 

Above 6M , both MFs differ significantly. The core MF features a power law with 
a flat slope, T = —0.26, indicating a very efficient production of high-mass stars in the 
Arches core. The sharp transition towards a flat MF above 10 M Q lets us speculate that a 
distinct physical process influences the formation of high-mass stars, argueing for primordial 
segregation. Such a scenario is consistent with the requirement of an exceptionally dense 
environment to form the most massive stars. The MF slope in the 2 nd annulus steepens to a 
Salpeter power law with T = —1.21 for M > 16 M Q , while intermediate mass stars are still 
overabundant as indicated by a flattened MF with T = —0.7 for 6 < M < 16 M Q . Before 
field contamination starts to dominate at 4M , the MF flattens to T ~ below 6M . 

A KS test comparing the incompleteness-corrected cumulative mass functions (CF) for 
M > 6 M yields a probability of less than 10~ 4 that the core and "field" (r > 13") CFs 
are drawn from the same distribution, and 0.03 of the core compared to the 2 nd annulus. 
The probability of a common origin of the 2 nd annulus and field CF is as high as 0.26, 
or 26%. The higher values derived for the 2 nd annulus indicate that the population at 
intermediate radii has some similarity with both the core and the field, but resembles the 
field much more closely. The KS test results support the core population to be distinct from 
mass distributions at all larger radii, and in view of the flat MF strengthens the high-mass 
bias/low-mass depletion in the Arches core. 

Indirect evidence for a low-mass truncated MF has been found in young, massive clus- 
ters in the center of the starburst galaxy M82 (Rieke et al. 1993). From the mass-to- 
light ratio, a low-mass cut-off at 2 — 3 M & is deduced assuming a normal high-mass IMF 
(Smith & Gallagher 2001). Below 3M Q , the low-mass MF in the Arches core is remarkably 
consistent with the area-scaled distribution of "field" stars in the cluster halo (r > 0.5 pc), 
possibly indicating the low-mass limit in the IMF in the GC. However, tidal stripping can as 
well be responsible for the re-distribution of low-mass stars into the cluster outskirts. Such 
an effect, if at work in the Arches cluster, could also be at work in M82. 
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5. Summary and Discussion 

We report direct evidence for a low-mass depleted MF in a near-by starburst cluster. 
The present-day MF of the Arches cluster displays a TO around 6 — 7 M & . A strong bias to 
high- and against low-/intermediate-mass stars is observed in the cluster core, which can be 
caused by rapid dynamical segregation or the selective formation of high-mass stars on the 
cost of intermediate mass members. Field stars appear to become dominant around a few 
solar masses, indicating that proper field star subtraction might cause a decreasing low-mass 
MF for M < 4 M & . The apparent decrease in the low- mass MF can either be caused by tidal 
stripping of stars below 6 M due to tidal disruption of the cluster in the GC tidal field, or 
a truncated IMF in the dense GC environment. 

Star formation theories predict that the characteristic mass depends on the initial con- 
ditions of the molecular cloud core and its environment. For a cloud of temperature T c \ 
with an embedded core T core , the critical mass for cloud collapse is the Bonnor-Ebert mass, 
M BE w 0.35(T core /10K) 2 • (n^T^/lO^rg/cm 3 )- 1 / 2 (Elmegreen 2000). In molecular clouds 
in the solar neighbourhood, typical values of nn 2 — 10 4 cm -3 , T c i = 100 K and T core = 10 
K result in a characteristic mass of ~ 0.3 M & . Star formation in the GC environment takes 
place under more extreme conditions. Increased temperatures and densities, turbulent pres- 
sure, and the influence of magnetic fields (Mouschovias 1991, Shu et al. 2004) may alter the 
outcome of the star formation process. In the Central Molecular Zone the average T core = 70 
K and cloud pressures are higher due to turbulence (Morris & Serabyn 1996). Assuming 
the thermal pressure is enhanced by a factor of 10 yields a Bonnor-Ebert mass of ~ 5M . 
Following Shu et al. 2004, the same characteristic mass arises when taking into account a 
magnetic field with B = 100/iG (Morris & Serabyn 1996). Hence the TO in the Arches MF 
can be understood as a natural result of the extreme star-forming environment in the GC. 

If a low-mass depleted IMF can be confirmed for the Arches cluster, this implies that 
a low- mass truncated and/or flattened IMF could be more adequate in intense star- forming 
environments, especially when using stellar synthesis models to derive the underlying mass 
distribution and star-formation history of star clusters and galaxies. While the integrated 
luminosity will still be determined by high-mass stars, the underlying population will con- 
tribute significantly less to the build-up of a Galactic field population of low-mass stars 
surviving throughout a Hubble time. Metal enrichment will appear enhanced when com- 
pared to the surviving low-mass population after a timescale of a few 100 Myr to a few Gyr, 
as the fraction of high-mass stars expected from a standard IMF will underestimate the true 
number of massive stars that formed and evolved in a galaxy. These implications have to be 
taken into account in the derivation of the star-forming histories and evolution of external 
galaxies, and especially in the extreme star-forming environments of distant galaxies in the 
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early universe. 

The question whether the Arches IMF was indeed truncated below a few solar masses, 
whether the GC environment caused an unusually high TO mass as compared to the low- 
mass TO at 0.2-0.5 M Q observed in the Trapezium (Muench et al. 2002; Hillenbrand 1997), 
or whether the cluster's dynamical evolution in the GC tidal field can entirely account for 
the present-day MF can only be solved by investigating a large area around the cluster and 
detailed field observations at comparably high spatial resolution as the NACO data presented 
here. 

AS acknowledges DGDF grant support from the ESO, and wishes to thank MPIA and 
Hans- Walter Rix for being an unusually supportive thesis host. We thank our referee, Simon 
Portegies Zwart, for suggestions that helped clarify and enhance the paper. 
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Fig. 1. — NACO combined HK image of the Arches cluster. The field of view is 24 x 24 
with a final resolution of 0.'084. North is up, East left. 
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Fig. 2. — Differential extinction corrected CMD of the Arches cluster. The 2 Myr Geneva 
isochrone used to derive the MF is shown as the solid central line. The color cut used to 
reject foreground and background contamination (enveloping lines) widens according to the 
increasing photometric uncertainty at fainter magnitudes. The 50% incompleteness limits 
in the cluster core and the 2 nd annulus (0.2 < r < 0.4 pc) are shown as dashed lines, and 
masses are labeled as derived from the isochrone displayed. The radial selection to r < 10" 
or 0.4 pc minimizes field star contamination. 
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Fig. 3. — Present-day MF as derived from a 2 Myr Geneva isochrone applied to the CMD 
shown in Fig. 2. A bin width of A log M/M = 0.2 and bin shift of 5 log M/M© = 0.02 
was used to obtain individual number counts. The gray enveloping lines denote the Poisson 
uncertainty for each point. The slopes and fitting uncertainties are derived from a weighted, 
linear least-squares fit between 6 < M < 60 M (dashed line). The high-mass MF is with 
T = —0.86 moderately flattened with respect to a Salpeter IMF (r = —1.35), indicating a 
bias to high-mass objects, and displays a TO at log M/M Q = 0.8, M = 6M Q . An upper 
limit to the possible field star contamination (dash-dotted line) derived from the outer NACO 
field, indicates that field stars are dominant below logm = 0.6 or 4M Q . 
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Fig. 4. — Mass segregation and turn-over in the Arches cluster core. The core MF (r < 0.2 
pc) (left) is heavily biased to high-mass stars, with a pronounced TO at logm = 0.84 or 
7 M & . The linear, least-squares fit (dashed line) shows that the MF steepens from T ~ —0.26 
(12 < M < 60 M Q ) to a Salpeter-like slope of T = -1.21 (16 < M < 60 M Q ) at the high- 
mass end and T = —0.69 (6 < M < 16 M ) for intermediate masses in the 2 nd annulus 
(0.2 < r < 0.4 pc), while the TO is still present around 6M . The area-scaled MF of 
the outer NACO field (r > 0.5 pc) indicates that the TO becomes more pronounced when 
field data are available. The similar TO mass in both annuli contradicts expectations from 
dynamical segregation and supports a low-mass depleted IMF in Arches. 



